Incommensurate nematic fluctuations in the two-dimensional Hubbard model 
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We analyze effective d-wave interactions in the two-dimensional extended Hubbard model at weak 
coupling and small to moderate doping. The interactions are computed from a renormalization 
group flow. Attractive d-wave interactions are generated via antiferromagnetic spin fluctuations in 
the pairing and charge channels. Above Van Hove fllling, the d-wave charge interaction is maximal 
at incommensurate diagonal wave vectors, corresponding to nematic fluctuations with a diagonal 
modulation. Below Van Hove fllling a modulation along the crystal axes can be favored. The nematic 
fluctuations are enhanced by the nearest-neighbor interaction in the extended Hubbard model, but 
they always remain smaller than the dominant antiferromagnetic, pairing, or charge density wave 
fluctuations. 



I. INTRODUCTION 

Numerous layered correlated electron compounds ex- 
hibit competing instabilities, which are driven by a vari- 
ety of dynamically generated effective interactions. This 
feature is already borne out by the simplest model of a 
single layer in such systems, the two-dimensional Hub- 
bard model. Near half-filling, a whole zoo of enhanced 
fluctuations is generated in that model, not only in the 
doped Mott insulator regime at strong coupling,^ but also 
for a weak bare interactioni^"— In addition to the pro- 
nounced antiferromagnetic and d-wave superconducting 
fluctuations, an attractive charge forward scattering in- 
teraction with a dx2 ^y2-wa,ve symmetry was found,^'^ If 
strong enough, such an interaction can trigger a d-wave 
Pomeranchuk instability leading to a nematic state with 
broken orientation symmetryi^ 

Nematic order with a d^^-y^ symmetry has been ob- 
served in several layered correlated electron compounds. 
A nematic phase with a sharp phase boundary has 
been established in a series of experiments on ultrapure 
Sr3Ru207 crystals in a strong magnetic fieldiiS Elec- 
tronic nematicity has also been observed in the high 
temperature superconductor YBa2Cu30y in transport 
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and neutron scattering 
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Nematic tendencies compete with other instabili- 
ties. In a functional renormalization group (fRG) 
analysis, with forward scattering and other interac- 
tion channels treated on equal footing, nematic fluc- 
tuations in the two-dimensional Hubbard model were 
found to be weaker than antiferromagnetism and d-wave 
superconductivity.— However, this does not necessarily 
prevent a nematic state, since nematic order may coexist 
with antiferromagnetism^^ and superconductivity i^^'^^ 
Moreover, the fRG calculations are reliable only in the 
weak coupling regime. A pronounced nematicity has 
been obtained by dynamical cluster mean-field calcula- 
tions for the two-dimensional Hubbard model at strong 
coupling ji^ 

In a recent analysis of secondary instabilities generated 
by antiferromagnetic fiuctuations in a two-dimensional 
electron system with a nearly half-filled band, Metlit- 
ski and Sachde v^^'^" found a tendency toward forma- 



tion of a modulated nematic state. In that state the ne- 
matic order oscillates across the crystal, with generally 
incommensurate wave vectors that point along the Bril- 
louin zone diagonal and connect antiferromagnetic hot 
spots with collinear Fermi velocities. At the same wave 
vectors the bare d-wave polarization function of tight- 
binding electrons on a square lattice has pronounced 
peaks, which also indicates that a modulated nematic 
state is favored over a homogeneous one.^^ Upon reduc- 
ing the electron concentration the hot spots move toward 
the saddle points of the dispersion and the modulation 
vector shrinks, until it vanishes at Van Hove filling. Be- 
low Van Hove filling a modulation with a wave vector 
along the crystal axes can be favored. ^^ 

In this article we analyze the tendency toward forma- 
tion of a modulated nematic state in the two-dimensional 
Hubbard model and its extension with a nearest-neighbor 
interaction. To this end we compute the effective two- 
particle interaction from a fRG flow equation and com- 
pare the nematic channel to antiferromagnetic, charge 
density, and pairing interactions. All channels are treated 
on equal footing and the approximations are controlled at 
weak coupling. We find that away from Van Hove filling a 
modulated nematic instability can indeed be favored over 
a homogeneous one, but remains in any case subleading 
compared to antiferromagnetism and d-wave supercon- 
ductivity. The nearest-neighbor interaction present in 
the extended Hubbard model enhances the nematic fluc- 
tuations. 

The paper is structured as follows. In Sec. II we intro- 
duce the extended Hubbard model and some notation. 
In Sec. HI we describe the fRG flow equation and specify 
the parametrization of the effective two-particle interac- 
tion. Results for the effective interaction are presented 
and discussed in Sec. IV. A summary with the main con- 
clusions follows in Sec. V. 



II. EXTENDED HUBBARD MODEL 

We will analyze the extended Hubbard model on a 
square lattice with a nearest neighbor hopping amplitude 
t > and a next-to-nearest neighbor hopping t' < 0. In 



standard second quantization notation, the correspond- 
ing Hamiltonian reads 



H =y e^al^ak^ 
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with a tight-binding dispersion of the form ek = 
— 2i(cos fc^j-t-cos ky)—At' cos k^ cos ky. The Hubbard inter- 
action C/ > is an on-site repulsion of electrons with op- 
posite spin. The operator njo- counts electrons with spin 
orientation a on the lattice site j , and jij = Uj^ + rij^. 
The second interaction term with y > is a density- 
density repulsion between nearest neighbor lattice sites. 
The lattice sum over nearest neighbors {j,j') is defined 
such that each bond is taken into account only once. The 
nearest neighbor interaction favors charge order. In mean 
field theory there is a first order transition between spin 
and charge density wave order as a function of V^^^ 
We shift the chemical potential /i by 4t' such that for 
^ = the Fermi surface contains the saddle points (tt, 0) 
and (0, tt), where the gradient of ek vanishes. The corre- 
sponding filling is referred to as Van Hove filling. 

In a functional integral formalism, a one-band model 
of interacting fermions is described by a bare action of 
the formal 
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dfc }^ ^Ak)[-iko + (ek - ^J-)]Mk) + Vo[*] 



(2) 
with anticommuting Grassmann fields tp and ip. In ad- 
dition to momentum k and spin a the fields V'cr(^) ^nd 
ipa{k) also depend on Matsubara frequencies /cq. The 
variable k — (fco,k) collects frequencies and momenta, 
and / d/c is a short-hand notation for J ^- J -^^ ■ We 
consider only the case of temperature zero in this paper, 
so that the Matsubara sums are integrals. The symbol ^ 
denotes a dependence on both ip and ip. For the extended 
Hubbard model, the bare two-fcrmion interaction can be 
written as 

"^o[*] = 2 / dpidp2dp3 V'o(pi,P2,P3) 

X YfPa{PlHcT'iP2)lpa'{P3)lpa{P4), (3) 
fT,cr' 

where Pi = pi + P2 — Ps, is fixed by momentum and fre- 
quency conservation, and 

^o(Pi,P2,P3) = C/ + V^.9(P2 -Pa), (4) 

with g{q) = 2(cos(7j; -I- cos^y). 



a smooth frequency cutoff A, corresponding to a regular- 
ized bare propagator 
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The effective interaction V"^ on scale A interpolates 
smoothly between the bare interaction Vo at A = oo 
and the final effective two-particle interaction obtained 
in principle for A ^> 0. The exact flow of V^ is deter- 
mined by a differential flow equation, which involves ef- 
fective ?n-particle interactions of arbitrary order m^ We 
use a weak-coupling truncation of the exact flow equa- 
tion, where ?7i-particle contributions with m > 3 and 
self-energy feedback to the flow of V"^ are neglected. Cor- 
rections to this approximation arc of order (V^)'^. The 
truncated flow equation has the schematic form 



d 



A 



1 

2*^ 






d^v' 






d'^v' 



9*2 U 9^2 



(6) 



where Gq 



d /~<A 
dA'^0- 



Graphically this corresponds to two 
vertices connected by two regularized propagators (one 
differentiated) forming particle-particle or particle-hole 
bubbles. 

We neglect the frequency dependence of the effective 
interaction, which is irrelevant in power-counting at weak 
coupling. To obtain an efficient parametrization of the 
momentum dependence we make an ansatz^^ 



V\^ ^ Vo[*] + V,^[*] + V,t[*] + V^[*] 



(7) 



which is a decomposition into spin, charge, and pairing 
channels. The spin or magnetic channel describes the 
interaction of spin operators 
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*] = -E / dg A/„^q)5]5i'^)(9)5,i'^H-g) (8) 
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Here r*^"-' are Pauli matrices, and /i(k) = 1 and /2(k) — 
cos fc^,— cos ky are s-wave and d-wave form factors, respec- 
tively. The coupling functions M^(q) evolve in the RG 
flow. If Af;^(q) develops a peak at (close to) q = (tTjTt) 
which diverges upon lowering A, the system tends to (in- 
commensurate) antiferromagnetic order. 

A similar parametrization is chosen for the charge 
channel 



v^m 



1 ^ 
-Y 

n=l 



dq<(q)A^„(q)7V„(-g), (9) 



III. WEAK COUPLING FUNCTIONAL RG 

We compute the effective two-particle interaction of 
the extended Hubbard model from a functional RG flow.^ 
To this end we regularize the bare action by introducing 



where Nn{q) = j '^kY.^'iP^{k)%P„{k + g)/„(k + f ) are 
density operators. An enhancement of K^{c() at a fi- 
nite wave vector q indicates a tendency toward forma- 
tion of a charge density wave. In this work, we pay par- 
ticular attention to the coupling function Xj (q) near 
q = 0, which, if sufficiently large, can drive a d-wave 



Pomeranchuk instability leading to a nematic state. At 
Q = (tt, tt) the coupling i^2(Q) will not be very pro- 
nounced, because for k sitting on a Van Hove point, 
k + ^ points to a zone diagonal where the form fac- 
tor /2 is zero. Note that (iV2(Q)) with Q = (0,Q) dif- 
fers from the familiar d-density wave order parameter— 

The pairing channel describes the interaction 

V^m^~J2 /d<zi?^(q)X„(g)X„(q) (10) 

with spin singlet Cooper pairs 

X{q) = Jdk ^^{k)lp^{q - fc)/„(f - k) , 

x{q) ^ Jdk Mk)Mq - k)fn{% - k) . 

The form factors determine the symmetry of the gap. A 
diverging (upon lowering A) peak of -Df (q) at q = 
indicates an s-wave pairing instability, and a diverging 
peak of -02^ (q) at q = a d-wave pairing instability, so 
that a superconducting state is favored. Since both form 
factors are even functions of momentum we only consider 
singlets in the spin dependence. 

The restriction to s-wave and ti-wave form factors in 
our ansatz is motivated by the dominance of s-wave and 
d-wave interactions observed in all previous fRG studies 
of Hubbard-type models in the parameter range consid- 
ered here.— 

Substituting the ansatz Eq. ([7]) into the flow equation 
generates a variety of different one-loop graphs, not all 
being of the form of the ansatz. We will use a suitable 
approximate projection which was developed in Ref. |25| . 
First, the graphs are sorted according to the transfer mo- 
mentum that propagates through the fermion loop. As 
long as the coupling functions do not have pronounced 
peaks this is the main momentum dependence and the 
graphs can be assigned to one channel. Since the form 
factors /i and /2 are orthonormal functions on [— 7r,7r]^, 
the graphs in each channel can be projected to the flow 
of coupling functions by integration. The resulting flow 
equations can be found in Ref. 25 and have been shown 
to yield a reasonable approximation to the RG flow. 

The above truncation of the RG is a weak coupling ap- 
proximation. For high scales, U / K is a small parameter, 
which allows to truncate at one-loop order. Generically 
couplings start to grow when the scale is descreased. In 
an intermediate scale regime there are phase space argu- 
ments such that the one-loop flow remains valid although 
the couplings are not small anymore. ^^ Eventually, the 
coupling functions become too large for certain momenta 
and the truncation becomes unreliable. We then stop the 
flow and interpret the rapid growth of the leading cou- 
pling as an instability indicating a corresponding ordered 
state. Specifically, we stop the flow when the maximum 
of the coupling functions exceeds T^nax — 20t. The stop- 
ping scale A* is an upper bound for the critical scale, at 
which interactions diverge. 



IV. EFFECTIVE D-WAVE INTERACTIONS 

We now present results for the effective d-wave interac- 
tions in the charge and pairing channel in the presence of 
strong antiferromagnetic correlations. Throughout this 
section we choose a relatively weak Hubbard interaction 
U = it. 

The leading instabilities are signalled by diverging cou- 
pling functions in the corresponding channel. The diver- 
gence of the full two-particle vertex is well captured by 
these coupling functions. However, to compare effective 
interactions quantitatively, one has to sum the contribu- 
tions from all terms in the decomposition Eq. ([7]), in- 
cluding also the bare interaction Vq. Contributions from 
the bare interactions and other channels can be signifi- 
cant in channels where the effective interaction remains 
bounded. 

The complete effective interaction in the various chan- 
nels can be extracted from the two-particle vertex 
r^i(T2;<T3<T4(Pi:P2;P3,P4), which is the antisymmetrized 
kernel of V"^, that is, 

^^[*] = 4 51 /dpidp2dp3r^i.,;.3.,(pi,P2;P3,P4) 



(Ti,...,cr4 
X -0^1 (Pl)V'<T2 (P2)V'<T3 (P3)V^<T4 [Pi] : 



(11) 



where P4 = pi + P2 — P3- We denote the static limit 
of F'^ by F^^^^.^^^^(pi,p2;p3,P4). The effective singlet 
pairing interaction is then given by 



Tpp' =rn;;T(P'"P;"P'P)' 



(12) 



and the effective charge interaction by 

rpp'^(q) = I E rLv'.(P' p'; p' - q. p + q) . (i3) 



where q is the momentum transfer. In a renormalized 
mean- field theory of spin singlet superconductivity based 
on F , only the pairing interaction Fpp', contributes.— 
Similarly, a renormalized mean-field theory for charge or- 
der with a modulation vector q involves exclusively the 
charge interaction Fpp,(q). 

The s-wave and d-wave components of the interactions 
are extracted by the projections 

7n^'^ =J^.^. fnip)fM r-/ (14) 

and 

7n''(q) = / (^1^ /"(p + i)/«(p' - f ) r^p'(q) 

(15) 
with n = 1,2. Inserting the channel decomposition 
Eq. ([7]), the d-wave components can be written as 
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(16) 



and 
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i^) = -^-iK^i<l) + 



(17) 



where the dots refer to fluctuation contributions from 
other channels. Complete expressions for 7^^^'^ and 
7^'^(q) in terms of the bare interaction and the cou- 
pling functions M^, K^ and D^ are given in the Ap- 
pendix. Note that positive coupling functions yield at- 
tractive (negative) contributions to the corresponding ef- 
fective interactions. 



A. Van Hove filling 

We first consider the case of Van Hove filling for various 
choices of t' < 0. For \t'\ < 0.237i the RG flow runs into 
an antiferromagnetic instability, that is, Mi{q) with q 
close to (7r,7r) reaches Vmax first. For \t'\ < O.llt the 
maximum of M-^ * (q) is situated at Q = (7r,7r), while it 
deviates from (tt, tt) for larger \t'\. However, this shift 
of q depends sensitively on the choice of A*, and the 
maximal value of M-^ * (q) differs only very little from its 
value at Q. For \t'\ > 0.237t the dominant instability is 
d-wave pairing, that is, -D^(O) is the largest coupling at 
A*. The transition from antiferromagnetism to d-wave 
pairing upon increasing \t'\ was found already in the first 
fRG studies of the two-dimensional Hubbard modelj^i^ 

In Fig. 1 we show results for the effective d-wave in- 
teractions at the scale A, as a function of t' /t. The an- 
tiferromagnetic (AF) and d-wave pairing (dSC) regions 
are separated by a vertical dashed line. The stopping 
scale A* decreases monotonically from 0.194i at i' = to 
0.013i at t' = -0.28i for V^ = 0, and from 0.162i at i' = 
to 0.007t at t' = ~0.28t for V = 0.74t. At Van Hove fill- 
ing, the d-wave charge interaction 73 ' * (q) is peaked at 
q = 0. For the plain Hubbard model {V = 0) it is always 
weaker than the d-waye pairing interaction 72 ' * , as ob- 
served already in Ref. [ij. A finite V enhances the d-wave 
charge interaction, but it remains small compared to the 
dominant antiferromagnetic and pairing interactions at 
small and large |i'|, respectively. 

In Fig. 2 we plot the effective interactions 73 ' * (0) and 
7!"'^* as functions of V at t' = -0.15t. From ^ = to 
V = 0.758i the stopping scale decreases moderately from 
A* = 0.078t to A, = 0.053i, but it then increases rapidly 
to A* = 0.656t at V = t. For increasing V < 0.758i, 
in terms of absolute values, the d-wave interaction in 
the density channel increases while the d-wave pairing 
interaction decreases. For V > 0.758i the leading insta- 
bility changes from incommensurate antiferromagnetism 
to a charge density wave, that is, the coupling function 
K^{Q) reaches Vmax first and defines the stopping scale 
A,. The critical value V = 0.758i is close to the mean- 
field transition point2^ at V = 11/4 = 0.75i. The drop 

of -72 '"^'(O) and -7^'"^* for V > 0.758t is due to a 
rapidly increasing A» (for increasing V) in the charge 
density wave regime. 




FIG. 1: (Color online) Effective d-wave interactions at the 
scale A* in units of t as a function of t' , for U — 3t a.t Van 
Hove filling. The electron density ranges from n = 1 at t' = 
to n = 0.75 at t' /t — —0.28. The pairing interactions 
71"^' * are plotted as circles and charge interactions 7^' * (0) 
as squares. Solid lines correspond to V = and dashed lines 
to V = 0.74t. In both cases the leading instability changes 
from antiferromagnetism to d-wave pairing at t' = — 0.237f 
(corresponding to n = 0.79). 




FIG. 2: (Color online) Effective d-wave interactions 7!*^' * 
(circles) and 7^' * (0) (squares) in units of t as a function 
of the nearest-neighbor interaction V, for t' = — 0.15i and 
U — 3t at Van Hove filling (n — 0.88). The leading instability 
changes from incommensurate antiferromagnetism to charge 
density wave at V = 0.758t. 
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FIG. 3: (Color online) Effective d-wave interactions at the 
stopping scale A» of antiferromagnetism plotted as functions 
of jj, for t' = — 0.15t, U = 3t, and V = 0. The density varies 
from n = 0.94 at ^/t = 0.1 to n = 0.81 at n/t = -0.11. Pair- 
ing interactions 7!*^' * are plotted as circles connected by a 
solid line. The most attractive d-wave charge interactions 
maxq |7^' * (q)| are plotted as squares connected by a dashed 
line, while the points corresponding to 7J' * (0) are connected 
by a solid line. Also shown are momenta q where — 7J' * (q) 
becomes majdmal, with components Qx plotted as open dia- 
monds and Qy as filled diamonds, both connected by a dashed 
line. Momenta for which the bare d-wave particle-hole bubble 
^ph ' (q) is extremal are shown similarly with components qx 
plotted as plus signs and qy as diagonal crosses. Stars are 
superpositions of plus signs and crosses for q^ = qy 



B. Away from Van Hove filling 

In the following we fix the next-to-nearest neighbor 
hopping amplitude to t' — — 0.15t and vary the chemical 
potential /i in a limited range around Van Hove filling. 
For this choice of parameters the RG flow runs into an an- 
tiferromagnetic instability, that is, M(^{q) with q close 
to (7r,7r) reaches Vmax first. For /i > 0.03i the maxi- 
mum of M]^ * (q) is at Q = (tt, tt), and for /i < 0.03t four 
maxima are found at {tt,tt ± 6) and (tt ± (5, tt), with 6 
increasing upon lowering /i. The dominance of antiferro- 
magnetic correlations is not affected by a nearest neigh- 
bor interaction V as long as I^ < U/A. At sufficiently 
low /i (large hole doping), d-wave pairing becomes the 
leading instability.2r— Here we do not enter the d-wave 
superconducting region. 

Results for the effective d-wave interactions 72 ' * and 
72 ' * in the antiferromagnetic background are shown for 
the plain Hubbard model {V = 0) in Fig. 3. We distin- 



guish between 72 ' *(0) (solid line) and maxq I72' *(q)| 
(dashed line). Above Van Hove filling (/i > 0) there are 
four maxima of — 72 ' * (q) at q = (±(5i, ±(5i) with 61 de- 
creasing as we approach Van Hove filling. The positions 
of these maxima are plotted in Fig. 3 as well. They are 
very close to the positions of the maxima of the static 
regularized d-wave particle-hole bubble 



K.i^) - / dp G^ipo,p)G^ipo, P + q)/2(p + f )' 



(18) 

at A = A*, which are plotted in Fig. 3 for comparison. 
The small deviations are not significant, that is, they 
may be due to the limited momentum resolution in the 
numerics. 

A large attraction j^' *(q) ^^ ^ wave vector q y^ 
would signal a tendency to form a modulated nematic 
order with a modulation vector q. Such a tendency 
was found in a recent analysis of secondary instabili- 
ties generated by antiferromagnetic fiuctuations in a two- 
dimensional metal by Metlitski and Sachdev i^^i^° The 
modulation vectors q are fixed by the position of the an- 
tiferromagnetic hot spots, where the antiferromagnetic 
fiuctuations couple most strongly to the electronic ex- 
citations near the Fermi surface. Subsequently it was 
shown that the same modulation vectors q = (±(5i, ±(5i) 
emerge as the momenta where the static d-wave particle- 
hole bubble is maximal. ^^ For a tight-binding dispersion 
with hopping amplitudes t and t', the modulation is given 



Si = 2 arccos . / 1 



4|t' 



(19) 



Note that ^ > has been shifted by 4t' such that Van 
Hove filling corresponds to 11 — 0. The maxima of the d- 
wave interaction — 72 ' * (q) obtained from the fRG flow 
above Van Hove filling are also situated at diagonal mo- 
menta of the form q — (±(5i,±(5i), with 61 given by 
Eq. (I19p within the numerical resolution. In agreement 
with Metlitski and Sachde v^^'^" we find that in the an- 
tiferromagnetic regime the d-wave density interaction is 
of the same order of magnitude but smaller than the d- 
wave pairing interaction. However, while d-wave pairing 
competes with antiferromagnetism at larger doping, the 
d-wave density interaction remains rather small. 

Slightly below Van Hove filling the particle-hole bubble 
is very flat. In the numerical solution of the RG flow we 
find the maximum of —72' *(q) at q = for — 0.08i < 
/i < 0. For lower electron filling, ^ < — O.OSi, there 
are again four degenerate maxima, but now at momenta 
along the crystal axes (0, ±^2) and (±^2,0). Overall, the 
strength of the effective d-wave interaction in the density 
channel is weaker than in the d-wave pairing channel and 
is maximal around Van Hove filling. 

A nearest neighbor repulsion is expected to enhance 
the d-wave charge correlations. Indeed, in a mean-field 
study restricted to spatially homogeneous solutions it was 
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FIG. 4: (Color online) Effective d-wave interactions and mod- 
ulation vectors as in Fig. 3, but for a finite nearest neighbor 
interaction V = 0.7 At. 



found that the nearest neighbor interaction in the ex- 
tended Hubbard model can generate a nematic statei^ 
In Fig. 4 we show results for the effective d-wave in- 
teractions for the extended Hubbard model with a finite 
nearest-neighbor interaction V ~ 0.74i. The other pa- 
rameters are the same as in Fig. 3. The d-wave den- 
sity interaction is enhanced by the presence of V, while 
the pairing interaction is reduced compared to the case 
V = 0. The favorite nematic modulation vectors are not 
affected significantly by V, but the peaks of 73' *(q) 
at those wave vectors are more pronounced in Fig. 4 
compared to the case V — plotted in Fig. 3. The 
stopping scale A^, decreases monotonically from O.lObt 
at n/t = 0.1 to 0.036i at fi/t = -0.11 ioi V = (Fig. 3), 
and from 0.076i at ^/i = 0.1 to O.Q19i at fi/t = -0.11 
for V = 0.74t (Fig. 4). 

We focused on the case of small or moderate dop- 
ing. Further away from half-filling, modulated ferro- 
magnetic fluctuations with incommensurate wave vectors 
have been found in mean-field and functional RG studies 
of the two-dimensional Hubbard model.—"— 



C. Mechanism 

To gain some analytical understanding of the above 
results, we now discuss the structure of the flow equations 
for the effective d-wave couplings. The flow equation for 
the d-wave coupling in the density channel can be written 



j25 

d 
dK 



K^i^) 



dA 



[K^{<l) + V 



a 



K\ 



(20) 



where a^ is a contribution from s-wave channels given 

by 



A 



d^k 



(cOSfca; -t- COS ky) 



[-2Z?f (k) + 3Afi^(k) + i^f (k)] . (21) 



The function i^j>(q) originates from d-wavc couplings 
D2, M2, and i^2 ■ We do not write this function here,— 
but it was taken into account in our numerical solution 
of the flow equations. As long as the d-wave interac- 
tions are small, contributions from F^{q) are negligible. 
Then the only momentum dependence is generated by 
a scale derivative of the regularized d-wave particle-hole 
bubble Eq. (|18p . The initial condition at a high scale Aq, 
derived from second order perturbation theory, is also 
given by this bubble K^°{q) = -V'^<^l;!"'{q) > 0. This 
makes plausible why the maxima of ii'^(q) and — 72 ' (q) 
are given by the extrema of the regularized particle-hole 
bubble. 

Even if there is no d-wave coupling initially (for V — 
0), K2 is generated in the flow similarly to the mecha- 
nism that generates an attraction leading to d-wave pair- 
ing. Close to an antiferromagnetic instability the mag- 
netic coupling function obeys M^(Q) » Af^(O) with 
Q = (tt, tt). Via integration with the cosine this leads to 
a sizeable negative contribution in Eq. (PT|) . Then the 
square bracket in Eq. (j20p remains non-zero for K2 > 
and increases with the latter. Since ^A^ph (q) > and 
the scale A decreases, the coupling function A'^(q) builds 
up in the flow. It is not essential whether M^ (k) is max- 
imal at or close to (tt, tt). The modulation of the d-wave 
coupling is not tied to an incommensurability in the mag- 
netic interactions. 

In the same way AT^ is generated by the density chan- 
nel if A'^(Q) :^ A'^(O), that is, if the system is close to a 
charge density wave instability (CDW). The latter is the 
case for larger V, which also helps to generate AT^ via the 
constant V in Eq. d^U]) and via a larger initial condition. 
Therefore we find a larger effective d-wave interaction in 
the density channel if the nearest neighbor interaction V 
is increased, as in Figs. 2 and 4. 

The flow equation for D^ is given by2^ 



dA 



^^(q) 



dA 



(22) 



where $2,A(q) = /dpG^(po, p)G^(-Po,q-p)/2(f -p)' 
is the static regularized particle-particle bubble with d- 



wave form factors and 

1 /• d^k 

4 



a 



J '^ ^^°'^" + ^°'^^^ [3A^i''(k) - i^f (k)] 



(23) 



Again we neglect the function F^ (q) for simplicity, which 
is justified for small K2 and M^. The initial condition 
at a high scale Aq , obtained from second order perturba- 
tion theory, is given by £'^°(q) = V^^l^f^iq) > and is 
dominated by the constant —V in Eq. ((22)) . Because 
dA^l'J^{0) < and A decreases, the coupling D^iO) 
would saturate at V if a^ — 0. The d-wave pairing can 
only become large if a^ is negative. Like for K2 this 
is the case in an antiferromagnetic background, which is 
the main mechanism for d-wave superconductivity in the 
Hubbard model. However, because of the minus sign in 
front of K^Ck) in Eq. (P^ . the vicinity of charge density 
order is counterproductive for the evolution of d-wave 
pairing. 

It is instructive to compare the contribution from mag- 
netic interactions to the generation of d-wave couplings 
described above to the hot spot mechanism discovered 
by Metlitski and Sachdevi ^^i^" They consider the situa- 
tion at a quantum critical point on the phase boundary 
of a commensurate antiferromagnetic ground state. As a 
consequence, their spin fluctuation propagator is strongly 
peaked at (tTjTt), and electronic excitations around hot 
spots dominate. The nematic modulation vector is deter- 
mined by the distance between hot spots with coUinear 
Fermi velocities. Metlitski and Sachdev find a degen- 
eracy between the d-wave couplings in the pairing and 
charge channels, which is broken only by the Fermi sur- 
face curvature. This degeneracy is present also in our 
flow equations Eqs. (pij) and (P5|) . where M^^Ck) con- 
tributes exactly equally to aj^ and a^. The Fermi sur- 
face curvature lifts the degeneracy because it reduces the 
size of the particle-hole bubble at the nematic modula- 
tion vector compared to the particle-particle bubble at 
q = 0. Singular self-energy contributions at the quan- 
tum critical point suppress the curvature effects, so that 
the nematic correlations are more pronounced than away 
from criticalityJ^ 



V. CONCLUSION 

We computed effective d-wave interactions in the two- 
dimensional extended Hubbard model at small to moder- 
ate distance from half-filling by using a functional renor- 
malization group fiow. In addition to the well-known 
attraction in the d-wave pairing channel, an attractive 
interaction in the d-wave charge channel is generated, in 
agreement with early fRG calculations.^ If strong enough, 
that interaction could induce a d-wave Pomeranchuk in- 
stability leading to nematic order. However, in compar- 
ison to the dominating antiferromagnetism and d-wave 
pairing, the d-wave charge interaction is found to be rel- 
atively weak. It becomes most pronounced at Van Hove 



filling and can be enhanced by a nearest neighbor inter- 
action V. 

The d-wave charge attraction is not necessarily max- 
imal at q = 0. Above Van Hove filling, we find four 
degenerate peaks of 72' (q) at diagonal wave vectors 
q = (±(5i, ±(5i) connecting hot spots with collinear Fermi 
velocities, corresponding to the modulated nematic cor- 
relations discovered by Metlitski and Sachdev i^^i^° In 
the fRG calculation, the peaks are essentially determined 
by the structure of the d-wave particle-hole bubble, dis- 
cussed in detail in Ref. [2l|. Below but close to Van Hove 
filling 1 72' (q)| is quite flat with a shallow maximum at 
q = 0. Further decreasing /i one obtains maxima along 
the crystal axes, again corresponding to extrema in the 
d-wave bubble. 

In any case, the nematic fluctuations in the weakly 
interacting extended Hubbard model seem to be gener- 
ically small compared to the dominant channels, which, 
depending on parameters, are antiferromagnetic, d-wave 
pairing, or charge density fluctuations at small or mod- 
erate doping. The nematic tendency observed experi- 
mentally in cuprates is thus probably a strong coupling 
phenomenon j "'^^1 ^^ possibly associated with antiferromag- 
netic quantum criticalityJ^^ For a microscopic theory 
of the nematic order in Sr3Ru2 07 a multi-band model^ 
seems to be required. 



Acknowledgments 

We are grateful to T. Holder, A. Katanin, M. Metlitski, 
and H. Yamase for valuable discussions. 



Appendix A: Expressions for effective interactions 

In this appendix we present explicit expressions for the 
s-wave and d-wave components of the effective charge and 
pairing interactions in terms of the bare interaction and 
the coupling functions. 

The s-wave component of the effective interaction in 
the singlet pairing channel as defined in Eq. P^ can be 
written as 



SC,A TT 

7i ' = [/ 



1 ^ r d^k 



(Al) 



3AC(k)-X;:(k) 



Similarly, the d-wave component is given by 

72 

hiik)+h2ik) 



1 /• d^k r 

+ 2 



(A2) 



(27r)2 

X 3Af2^(k) - Af (k) 

with a^ from Eq. (P^ . The functions /i„ are defined as 
hi{k) = i(cosfca;-|-cosfcj,) and /i2(k) =cos^cos-^. The 



first line is the main contribution and its square enters with a^ from Eq. ((2T|) . The first three terms enter the 
the flow equation of the coupling function Z)^ in Eq. (P^. flow equation for K2, see Eq. (pil)) . The remaining terms 



In the charge channel we allow for a finite transfer 
momentum q in the effective interaction between density 
pairs. Its s-wave component reads 



lt'\^) = \u + V9{ci)-\Kt{ci) 



c,A 



(A3) 



are subleading and are given by 
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c,A 



qx 



^2 (q) = 02,1 ^i(q) + 02.2 cos — cos 
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(A6) 



where 



c,A 



•'l.n 



d^k 

(2^)2 



2i?^(k)+3A/„^(k)+<(k) 



(A4) 



The d-wave component is given by 

72^^q) = ~\v - ^i^2^q) + \<^K + ^2^^q) , (as) 



with 



c,A 
^2,n 



d^k 



/i„(k) 



(2^)2 
X - 2D^ (k) + 3 Afa"^ (k) + i^^ (k) 



(A7) 
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